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Abstract 


Drought is one of the most frequent natural disasters to occur in north-central 
Mexico and its identification and characterization are crucial for mitigating its 
effects. This study analyzed the evolution of meteorological drought over the 


period 1961-2012 in the state of Zacatecas, Mexico, using three multiscale 
indices: the Standardized Precipitation Index (SPI) based on precipitation 
data, the Reconnaissance Drought Index (RDI), and the Standardized 
Precipitation-Evapotranspiration Index (SPEI) based on precipitation and 
evapotranspiration data. The temporal pattern of the drought was analyzed at 
three time scales (3, 6 and 12 months). The indices were compared using the 
Pearson and Spearman correlation coefficients. The main dry periods were 
1982-1983, 1998-2000, and 2010-2012. At a 3-month scale, the SPI detected 
the lowest number of drought events (41), with an average duration of 2.5 
months; and at 6- and 12-month scales, the RDI detected the lowest number 
of drought events, 32 and 16, respectively, with an average duration of 3.0 
and 6.1 months, respectively. The SPEI detected the greatest number of 
extreme events at the 3-month scale, and the RDI and the SPI did so at the 
6- and 12-month scales, respectively. For the three scales of analysis, the SPEI 
had the greatest average frequency percentage of severe droughts. The SPI 
vs RDI had the highest correlation (> 0.92), and the SPI vs SPEI had lowest 
correlation. The areas that were most affected by droughts were the 
northwestern, central, and southern regions of Zacatecas, where the principal 
rainfed and irrigated agriculture areas of the state are located. The findings 
can primarily help decision-making regarding cropping patterns and the 
optimal allocation of water in Irrigation District 034 in the state of Zacatecas. 


Keywords: Meteorological drought, SPI, RDI, SPEI, semiarid regions, severity 
and duration of droughts, potential evapotranspiration, Pearson and Spearman 
correlation coefficients. 


Resumen 


La sequía es uno de los desastres naturales más frecuentes en el centro norte 
de México, y su identificación y caracterización son cruciales para tratar de 
mitigar sus efectos. Se analizó la evolución de las sequías meteorológicas en 
el estado de Zacatecas, México, para el periodo de 1961 a 2012 utilizando tres 
índices multiescalares: el Índice de Precipitación Estandarizado (SPI), basado 
en datos de precipitación; el Índice de Reconocimiento de Sequía (RDI), y el 
Índice de Precipitación-Evapotranspiración Estandarizado (SPEI), basados en 
datos de precipitación y evapotranspiración. El patrón temporal de las sequías 
se analizó en tres escalas de tiempo (3, 6 y 12 meses). Los índices se 
compararon empleando los coeficientes de correlación de Pearson y de 


Spearman. Los principales periodos secos fueron de 1982-1983, 1998-2000 y 
2010-2012. En la escala de tres meses, el SPI detectó el menor número de 
eventos de sequía (41), con una duración media de 2.5 meses; y en las escalas 
de 6 y 12 meses, el RDI detectó el menor número de eventos de sequía, 32 y 
16, respectivamente, con duración media de 3.0 y 6.1 meses, 
respectivamente. En la escala de tres meses, el SPEI detectó el mayor número 
de eventos extremos, y en las escalas de 6 y 12 meses fueron el RDI y el SPI, 
respectivamente. En las tres escalas de análisis, el SPEI presentó el mayor 
porcentaje de frecuencia promedio de sequías severas. Los índices SPI vs. RDI 
presentan la correlación más alta (> 0.92), y el SPI vs. SPEI la correlación 
más baja. Las áreas más afectadas por las sequías son la región noroeste, 
centro y sur de Zacatecas, donde se ubican las principales zonas agrícolas de 
temporal y de riego del estado. Los resultados encontrados pueden ayudar 
principalmente en la toma de decisiones sobre los patrones de cultivo y la 
asignación óptima de agua en el Distrito de Riego 034 Estado de Zacatecas. 


Palabras clave: sequía meteorológica, SPI, RDI, SPEI, regiones semiáridas, 
severidad y duración de sequías, evapotranspiración potencial, coeficientes de 
correlación de Pearson y de Spearman. 
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Introduction 


Central and northern Mexico are located in the northern latitude high pressure 
belt. For this reason, these regions are arid and semi-arid, coinciding in 
latitude with the great deserts of Africa and Asia. Consequently, these regions 
do not have access to water in sufficient quantities and have historically been 
affected by droughts. In the 20th century, Mexico experienced four long 
periods of drought: 1948-1954, 1960-1964, 1970-1978, and 1993-1996, as 


well as a severe drought in 1998, which mainly affected the northern states 
(Cenapred, 2007). According to Escalante and Reyes (1998), the state of 
Zacatecas, located in north-central Mexico, is among the states whose 
agriculture and livestock was most damaged by the droughts that occurred 
between 1988 and 1995. 


Recently, several droughts of considerable magnitude have been registered, 
in the years 2000-2003, 2009, and 2011-2012. Between 2000 and 2003, 18 
of the 32 Mexican states suffered drought, and the northern states, including 
Zacatecas, were the most affected. Economic losses reached more than 1 800 
million pesos (188 million dollars). Nearly a million hectares of crops were 
affected, and more than 13 000 head of cattle were lost (Cenapred, 2002; 
Cenapred, 2003; Cenapred, 2004). In 2009, Mexico suffered the second worst 
drought in 60 years with damages reaching 3 081 million pesos (229 million 
dollars) and 384 540 hectares of crops and grasslands were impacted. 
Zacatecas was among the five states with the largest affected area (Semarnat, 
2012). In late 2011 and up to the middle of 2012, more than 60% of Mexico 
suffered the worst drought in the last 70 years, and its intensity was classified 
by the North American Drought Monitor (NADM) as severe to exceptional 
(Semarnat, 2012). Moreover, according to the United States Department of 
Agriculture (USDA), this drought caused economic losses to the agricultural 
sector of more than 16 000 million pesos (1 300 million dollars) (Herron, 
2013). The state of Zacatecas was the hardest hit by this natural disaster. 
Thus, it is important to have precise regional knowledge of droughts in the 
state of Zacatecas in order to establish mitigation measures. 


Droughts can have serious, lasting repercussions on human and natural 
systems, such as humanitarian disasters, economic losses, and stressed 
natural ecosystems (Touma, Ashfaq, Nayak, Kao, € Diffenbaugh, 2015). 
Droughts are driven by adverse climate variations and hydrological conditions 
(Li, Liang, Yu, € Acharya, 2014). Although they are a consequence of a natural 
reduction in the amount of rainfall in a period, their severity is evaluated in 
terms of their duration, the season in which they occur, the size of the affected 
area, and particularly, their impact on human and agricultural activities and 
the environment (Caparrini € Manzella, 2009). Droughts have been classified 
into four categories: meteorological, hydrological, agricultural, and 
socioeconomic (Samaniego, Kumar, € Zink, 2013). 


The complexity of defining a drought is due to the difficulty with quantifying 
its magnitude or severity, since it is typically identified based on its effects on 


different systems (water resources, agriculture, ecology, economic losses, 
etc.). However, there is no physical variable that can specifically measure the 
severity of a drought. The quantification of a drought generally uses drought 
indices (Asadi-Zarch, Sivikumar, 8 Sharma, 2015), which are based on data 
about precipitation, temperature, evapotranspiration, river flows, or other 
measurable hydrometeorological variables that quantify drought risk. Today, 
there are many indices for drought quantification. The most used worldwide 
are the Standardized Precipitation Index (SPI) developed by McKee, Doesken, 
8: Kleist (1993) and the Palmer Drought Severity Index (PDSI) developed by 
Palmer (1965). The SPI is a probabilistic approach to precipitation, while the 
PDSI analyzes soil water balance. Recently, the World Meteorological 
Organization (WMO) presented SPI as the universal drought index (Hayes, 
Svoboda, Wall, 8 Widhalm, 2011; WMO, 2012). The main advantage of the 
SPI is that calculations can be done at different time scales (usually 1, 3, 6, 9, 
12, 24, and 48 months), and hydrological variations of precipitation deficits 
can be represented. However, the SPI does have one disadvantage, namely, 
it uses only information on precipitation and does not consider other 
climatological variables that play an important role in the development of a 
drought event (Teuling et a/l., 2013). Recently, indices have been developed 
that incorporate evapotranspiration into the evaluation of droughts. For 
example, Tsakiris and Vangelis (2005) proposed the Reconnaissance Drought 
Index (RDT), which analyzes the ratio of accumulated precipitation to potential 
accumulated evapotranspiration (ETP) from a probabilistic standpoint. 
Vicente-Serrano, Begueria and López-Moreno (2010) developed the 
Standardized Precipitation-Evapotranspiration Index (SPEI), which consists of 
a monthly climate water balance. 


Because of the complexity of the drought phenomenon, none of the drought 
indices can be applied in all regions or to all natural systems. Therefore, to 
analyze droughts, more than one index should be used, in order to examine 
the sensitivity and precision of each index and to explore their behavior in 
specific conditions and for particular purposes. The characterization of a 
drought event involves identification of its start and end dates, duration, 
magnitude (severity), and intensity (McKee et al., 1993; Kavalieratou, 
Karpouzos, €: Babajimopoulos, 2012). Other characteristics of droughts include 
spatial coverage (area) and frequency of occurrence (Kirono, Kent, Hennessy, 
8: Mpelasoka, 2011) 


The objectives of this study were to (a) investigate the main characteristics of 
meteorological droughts in the state of Zacatecas during the period from 1961 
to 2012 using the multiscale drought indices SPI, RDI and SPElI; and (b) 
compare the temporal and spatial behavior of the drought indices used. 


Materials and methods 


Characteristics of the study area 


The state of Zacatecas is located between the coordinates 21% 01' 00” and 250 
07' 00” N and 1009 43" 00” and 1049 22' 00” W., Its area is 75 284 km?. The 
main climate types are arid and semiarid (73%) and temperate subhumid 
(17%). Mean annual temperature is 170C and average annual high and low 
temperatures are 30%C (May) and 3%C (January), respectively. The mean 
annual precipitation in Zacatecas is 510 mm, with 300 mm in the north and 
860 mm in the south. Most of the rainfall occurs in summer (75%), from June 
to September. The climate of Zacatecas limits agriculture, the main economic 
activity to which 1.74 million hectares is dedicated; 88% of this area is rainfed. 


Climatological information 


To identify the drought events, as well as to determine their duration, 
magnitude, intensities, and frequencies, this study used daily precipitation and 
temperature data series from 9 weather stations in the state of Zacatecas 


(Figure 1), taken from the CLICOM database (Conagua, 2014) for the period 
from January 1, 1961, to December 31, 2012. The geographic descriptors 
(coordinates and altitude) of the weather stations and the mean annual values 
corresponding to precipitation, P, and temperature, 7, are presented in Table 
1. 


Table 1. Weather stations and their corresponding mean annual precipitation 
and temperature. 


Latitu |Longit Precipita Temper 
de ude Altitud... tion 
Code |Name Period ature 
e (m) (0C) 
(N) (0) (mm) 
240 1020 1961- 
32052 San Rafael 35" 06' 2014 2012 262.8 14,4 
240 1030 1961- 
32028 Juan Aldama 16 23 1999 2012 433.0 18.0 
230 1020 1961- 
32001 Agua Nueva 47' 09' 1946 2012 349.1 17.7 
Sombrerete |230 1030 1961- 
2054 2 1.4 16.4 
pane (DGE) 38' 38' do 2012 39 E 
230 1030 1961- 
32018 El Sauz 16 06' 2096 2012 414.7 16.0 
220 1020 1961- 
32003 Calera 54 39" 2097 2012 426.5 15.6 
220 1030 1961- 
32030 La Florida 41 36' 1870 2012 590.3 16.4 
Guadalupe 220 1010 1971- 
32024 Victohia 23 49' 2132 2012 374.1 16.4 
ea 210 1030 1961- 
32032 La Villita 36' 20 1786 2012 782.3 20.4 
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Figure 1. Weather stations selected for drought analysis. 


Methodology 
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For this study, we used the three drought indices that are most commonly 
used and well-known worldwide: the Standardized Precipitation Index, the 
Reconnaissance Drought Index, and the  Standardized Precipitation- 
Evapotranspiration Index. Drought characterization was carried out in six 
stages: 1) selection and quality control of the climatological information; 2) 
selection of the scale of analysis; 3) calculation of potential evapotranspiration 
(ETP); 4) calculation of the SPI, RDI, and SPEI indices at time scales of 3, 6, 
and 12 months; 5) analysis of drought characteristics; and 6) comparison of 
the meteorological drought indices. Figure 2 schematically presents the 
methodological development of the study. 
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Selección de la información 


Pruebas de homogeneidad 
(t Student, Cramer, SNHT y Buishand) 


Complementación de datos faltantes 


(método de la relación normal) 


Etapa 2 Selección de la escala de análisis k 
(3, 6 y 12 meses) 
EApala Cálculo de la evapotranspiración 
(método de Thornthwaite) 
Etapa 4 Conformar las series acumuladas de: 


P (SPD o. (RDI) 


Ajustar las series D a una fdp 
(fdp Log-logística) 


D (SPEl) 


Ajustar las series P y QU a una fdp 
(fdp Gamma de dos parámetros) 


Características de la sequía 
(duración, magnitud e intensidad) 


Figure 2. Schematic diagram of the methodological development of drought 
analysis. 
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Selection and quality control of the climatological 
information 


A key stage in the analysis of hydroclimatological time series, as in the case 
of drought analysis, is the data selection and quality control. Selection of the 
climatological information used is primarily based on four criteria: (1) that the 
weather station is active, to be able to evaluate the droughts during the last 
few years; (2) that each weather station has precipitation and temperature 
records covering at least 40 years; (3) that the percentage of missing data is 
less than or equal to 15% in order to maintain data representativeness; and 
(4) that the stations are distributed throughout the study zone so as to have 
good spatial representation of the drought. 


For drought analysis, a preliminary quality control of the precipitation and 
temperature data series was carried out, checking for homogeneity of the 
climatological data series and filling missing data. To check for homogeneity 
of precipitation and temperature data, parametric and non-parametric 
statistical tests were applied to the annual series. The tests selected were the 
Student's t-test, the Cramer test (Mirza, 1997), the Standard Normal 
Homogeneity Test (SNHT), and the Buishand test (Hánsel, Mederios, 
Matschullat, Petta, 8 Mendoca-Silva, 2016). 


To select the method for filling in missing data, several methods were tested: 
polynomial interpolation, arithmetic average, inverse distance, and normal 
ratio (Tabios III € Salas, 1985; Dingman, 2002). The normal ratio had the 
best results (Equation (1)): 
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where fm is the missing datum to be estimated for the day, month, or year in 
the station of analysis m; g is the neighboring stations; n, is the number of 
neighboring stations considered in the analysis; p,, and p, are the annual mean 
precipitation at the station m with missing information and at the neighboring 
stations g, respectively; and p, is the precipitation registered at the 
neighboring stations on the day, month, or year in which the datum at station 
m was missing. This same method was also applied to the temperature data. 


Selection of the scales of drought analysis (k) 


The time scale, k, plays an important role in the study of drought 
characteristics. Short time scales are relevant because they are representative 
of meteorological and agricultural droughts. A time scale of 3 months gives a 
seasonal estimation of precipitation that reflects the moisture conditions on 
the short- and medium-term. A 6-month scale estimates precipitation on a 
medium-term scale and can be associated with reservoir levels and anomalous 
flows, while a time scale longer than 12 months reflects long-term patterns in 
precipitation that are likely related to flows, reservoir levels, and even 
groundwater levels, and permits the identification of hydrological drought. 


This study selected scales of 3, 6, and 12 months for the SPI, RDI, and SPEI 
indices (SPI-3, SPI-6, SPI-12, RDI-3, etc.), with the aim of obtaining short-, 
medium-, and relatively long-term perspectives of drought characteristics. 


Calculation of potential evapotranspiration (ETP) 
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To calculate the RDI and SPEI indices, potential evapotranspiration was 
estimated using the Thornthwaite (1948) method, as proposed by Tsakiris and 
Vangelis (2005) and Vicente-Serrano et a/. (2010). This method of calculating 
ETP (Equation (2)) is one of the simplest, since it only requires mean monthly 
precipitation data and the latitude of the site where ETP is to be estimated. 


mn (2) 


where ETP is the monthly potential evapotranspiration (mm/month); T is the 
monthly mean temperature (9C), and I/ is the heat index (Equation (3)), which 
is calculated as the sum of 12 monthly index values. 


0 o 


m is a coefficient that depends on the heat index 1. 
m=6.75x 107718 - 7,71 x 101? + 1.79 x 107?1+0.492 (4) 
K is a calculated correction coefficient as a function of latitude and month: 
TES o 
12 30 


where NDM is the number of days of the month and N ¡is the maximum number 
of hours of sunlight (Equation (6)). 
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Ni E US (6) 


and w, is the hourly sunrise angle (Equation (7)). 
w, = arccos(—tangtanó) (7) 


In which q is the latitude (radians) of the weather station and $ is the solar 
declination (radians) calculated as: 


. (21) 
$ = 0.4093 sin (E - 1.405) (8) 


where J is the average Julian day of the month. 


The Standardized Precipitation Index (SPI) 


The SPI uses precipitation as the only variable for the analysis, and estimates 
whether, in a given region and period, there is a deficit or excess of 
precipitation relative to normal conditions (Hayes, Svoboda, Wilhite, € 
Vanyarkho, 1999). The calculation of the SPI is based on long-term 
precipitation records. Guttman (1994) points out that preferably 50 to 60 years 
or more of records are needed, although other researchers (Wu, Hayes, 
Wilhite, £ Svoboda, 2005; Li et a/., 2014) recommend at least 30 years. The 
SPI is calculated by fitting an appropriate probability distribution function 
(pdf), generally the gamma distribution function, to cumulative monthly 
precipitation data for each time scale, k, and site of interest. The SPI values 
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are obtained by transforming the gamma distribution value into normal 
standard values. The gamma probability density function is defined by: 


aeo= cs TA forx>0 (9) 


1 
O 


where a is the shape parameter; b is the scale parameter; x is the cumulative 
monthly precipitation (x > 0), in mm, and TIí(a) is the gamma function. 
Parameters a and b can be estimated using several methods, and their values 
by the maximum likelihood method are (Haan, 1977): 


1 4A se 
E A (10) 


ES (11) 
a 


A= 100 -3Y x, (12) 


where n' is the number of precipitation data x, different from zero and x is the 
mean value of the x;, data. 


Integrating the gamma density function (Equation (9)), relative to x, we obtain 
the cumulative distribution function G(x) 


E E 
G(x) = Pr)” /b dx (13) 
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Considering t = x/b, Equation (13) is then reduced to the incomplete gamma 
function. 


1 x 
G(x) = ña) tLe-tdt (14) 


Given that the gamma function is not defined for x = 0, and the precipitation 
series can have zero values, then the cumulative probability can be calculated 
as (Angelidis, Maris, Kotsovinos, € Hrissanthou, 2012): 


H(x) = q + (1 q)G(x) (15) 


where H(x) is a function of the cumulative probability; q is the probability of 
having precipitation values equal to zero; and G(x) is the cumulative 
probability of the incomplete gamma function. 


If m is the number of zeros in the precipitation series, then q = m/n, where n 
is the total number of data in the precipitation series, according to the scale 
of analysis (k). When the series does not have precipitation values equal to 
zero, then q =0, and thus H(x) = G(x), and when x = 0, then H(0) = q. 


Finally, to define the SPI value, the cumulative probability function, H(x), is 
transformed into the standard normal variable, Z, which has a mean of zero 
and a variance of one (Edwards €: McKee, 1997), using the approach proposed 
by Zelen and Severo (1965): 


PEA , cc ttat us col 1 é 
Ae me ETRETEDET A AN dl (16) por (16) 


<H(x) < 0.5 
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Z =SPI = +(: 


< H(x) < 1.0 


Cy + Ct +c,t? nn 1 
14d t+d 24d e) (1-H(0) 


) for 0.5 (17) 


where c, = 2.515517, c, = 0.802853, c, =0.010328, d, =1.432788, d, = 0.189269, 
and d, = 0.001308. 


A drought event occurs when, at any time scale, SPI is continuously negative 
and reaches a value of -1.0 or less. It ends when the SPI becomes positive 
(McKee et al., 1993; Vrochidou € Tsanis, 2012). The SPI values can be easily 
compared simultaneously on the spatial and temporal dimensions (Lopez- 
Bustins, Pascual, Pla, 4 Retana, 2013). Based on the range of SPI values, 
drought events are classified from mild to extreme (Table 2). 


Table 2. Classification of SPI (McKee et al., 1993). 


SPI value Classification 

> 2.00 Extremely wet 

1.50 a 1.99 |Severely wet 

1.00 a 1.49 ¡MModerately wet 

0a0.99 Mild wet (close to normal) 
Mild drought (close to 

O a -0.99 normal) 

-1.00 a - 

1.49 Moderate drought 

-1.50 a - 

1.99 Severe drought 

< -2.00 Extreme drought 
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The Reconnaissance Drought Index (RDI) 


The RDI can be calculated at different time scales and its interpretation is 
similar to that of the SPI (Table 2) (Vicente-Serrano, Van der Schrier, 
Beguería, Azorin-Molina, € López-Moreno, 2015; Xu et a/., 2015). One of the 
most notable theoretical limitations of this drought index is that it is not valid 
when the ETP value is equal to zero, a very common situation in cold regions 
during winter (Vicente-Serrano et al., 2015). RDI is expressed in three ways: 
the initial value («,), normalized RDI (RDIn), and the standardized RDI 
(RDIst). The initial value (a, ) is presented in aggregated form using a monthly 
time scale and can be calculated on a monthly, seasonal, or annual basis 
(Asadi-Zarch, Malekinezhad, Mobin, Dastorani, € Kousari, 2011). Its values 
are calculated with the following equation: 


Di=1 Pi 


(0) 
Ó = 
de y ETP; 


SiS Ne (18) 


where P;, and ETP,, are the precipitation and potential evapotranspiration for 
month j of year i, respectively, in mm; N, is the number of years of available 
information; and k is the scale of analysis at which the RDI is calculated. The 
normalized and standardized forms of the RDI are (Tsakiris, Pangalou, € 
Vangelis, 2007): 


o _ 5 
RDI? = 5-1 (19) 
S, 
| M_5 
ROJO, =% Y (20) 
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where 8” is the arithmetic mean of 8? and y” are the natural logarithms of 
$; Whose mean and standard deviation are y, and 6,;, respectively. 


Different studies with data from different places and different time scales have 
found that, in all cases, the values of 6, followed either a gamma or a 
lognormal distribution, although in most of the cases the gamma distribution 
was more successful (Kousari et a/., 2014; Cai, Zhang, Yao, € Chen, 2015). 


Tigkas (2008) and Asadi-Zarch et al. (2011) have shown that the RDIst can 
be better calculated by fitting the gamma distribution to the distribution of 
$, frequencies, since this method tends to solve the problem of calculating 
RDIst for short time scales, which can include values of zero precipitation (9, = 
0), for which Equation (20) is not defined. 


In our study, the RDI was obtained by fitting the gamma distribution to 6, 
values, in a way similar to SPI, using Equations (10) through (17), and 
substituting the values of 6, instead of x. 


The Standardized Precipitation-Evapotranspiration Index 


The main limitation of the SPI is that it is based only on precipitation and 
ignores other variables that affect atmospheric water demand, such as 
temperature, wind speed, solar radiation, and vapor pressure deficit (McEvoy, 
Huntington, Abatzoglou, € Edwards, 2012). To overcome this limitation, 
Vicente-Serrano et al. (2010) developed the SPEI, which combines the PDST's 
sensitivity to changes in the evaporation demand, which is caused by 
temperature fluctuations and tendencies, with the SPI's simplicity of 
calculation and multi-scale nature (Vicente-Serrano et a/., 2010; Banimahd € 
Khalili, 2013; Li et a/., 2014). 


The SPEI is based on the monthly climatic water balance (precipitation minus 
potential evapotranspiration, called series D), which is calculated at different 
time scales and is mathematically similar to the SPI, but includes the role of 
potential evapotranspiration (López-Moreno et al., 2013). Its calculation 
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follows an approach similar to the one used to calculate the SPI, but with a 
three-parameter log-logistic probability distribution instead of the two- 
parameter gamma distribution (Figure 2). The SPI drought classification can 
also be used to evaluate the SPEI (Table 2). 


As a first step to calculate the SPEI, the difference between precipitation and 
potential evapotranspiration for month j of year i (series DY) is calculated, in 
mm. The series DY is then added in the time scale (k) of interest to obtain the 
series Dy?, in millimeters. 


DU = P;, — ETPy (21) 


where P;,; and ETP,; have the same meaning as in Equation (18). 


To calculate the SPEI, a three-parameter probability distribution needs to be 
used, since for a two-parameter distribution, the variable x has a lower limit 
of zero (0 < x < 0), whereas for the three-parameter distribution, x can take 
values in the range (y<x<o0), where y is the origin parameter of the 
distribution. For this reason, Vicente-Serrano et a/. (2010) proposed using the 
log-logistic distribution, since the D, series can have negative values. Its 
probability density function is: 


ta] ca 


where a, f, and y are the shape, scale, and origin parameters, respectively, 
for the D, values in the range of (y > D, < 00). 


The log-logistic distribution parameters can be obtained using different 
methods, of which the L-moments method is the most robust and the easiest 
to apply (Ahmad, Sinclair, € Werritti, 1988; Singh, Guo, € Yu, 1993). 
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where T(1 + 1/8) is the gamma function of (1 + 1/6) and w, is the moment of 
high probability of order s (s = 0,1,2): 


n 


an a 


i=1 


where x, is the ordered random sample (x, < x»,--*,< x,) of the D, values and n 
is the sample size. 


The probability distribution of the D, series, according to the log-logistic 
distribution, is given by (Vicente-Serrano et a/., 2010): 


HO)= ¡ de ST (97) 


Finally, to obtain the SPEI values, F(x) values are transformed into the 
standard normal variable, Z, substituting F(x) for H(x) in Equations (16) and 
(17). 


Comparison of drought indices 
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The SPI, RDI, and SPEI indices follow similar methodologies and their values 
have the same statistical significance; therefore, they are comparable. To 
compare the drought indices, we used the Pearson (r) and Spearman (r;) 
coefficients. 


Results and discussion 


Here, a drought event was defined as an event for which the values of SPI, 
RDI, or SPEI are equal to or less than -1.0. We determined the number of 
drought events and their duration and magnitude for scales of 3, 6, and 12 
months at each of the nine stations analyzed. 


A general analysis of the drought indices shows that, both locally and 
regionally, these indices follow the same temporal pattern for the period from 
1961 to 2012. According to the SPI, RDI, and SPElI indices, the wettest periods 
were 1965-1968, 1990-1992, and 2003-2005. The main drought events 
occurred in 1969-1970, 1982-1983, 1998-2000, and 2010-2012. 


The three indices analyzed showed consistent results in detecting the major 
drought events, in terms of their magnitude and duration (Table 3). The SPE, 
RDI, and SPEI series were very similar regardless of the time scale. These 
results indicate that for climate conditions with little interannual variability in 
temperature, the drought indices analyzed mainly respond to variability in 
precipitation. 


Table 3. Maximum values of drought characteristics for all stations and time 


scales. 
Index ¡Station | Duration Severity (Magnitude) Intensity 
Value Period Value Period Value Period 


22 


32052 16 on -8.426 | 12/2010 - 05/2011 | -1.483 bateo 
32028 (7 a -9.005 | 12/1988 - 06/1989 | -1.361 dos 
spPI-3  |32018 6 poa -9.696 | 12/2010 - 05/2011 | -1.713 pd 
32024 |6 oo -9.396 | 12/2000 - 05/2001 | -1.574 pd 
32032 (6 0 -8.936 | 12/2010 - 04/2011 | -1.787 Pa 
32052 (7 o -9.711 | 12/2010 - 06/2011 | -1.483 e 
32028 (7 os -9.109 |12/1988 - 06/1989 | -1.361 is 
RDI-3 |32018 |6 a -9.877 |12/2010 - 05/2011 | -1.713 pte 
32024 16 001 -9.414 | 12/2000 - 05/2001 | -1.574 a 
32032 16 RO -10.102 | 12/2010 - 05/2011 | -1.787 td 
32052 6 lan -9.775 [05/1974 - 10/1974 | -1.900 Dee 
32028 |5 pos -8.923 | 05/1980 - 08/1980 | -2.231 Deia 
SPEI-3 [32018 |5 ape -7.463 | 05/1962 - 09/1962 | -2.048 De tico 
32024 |4 re -7.247 | 04/1989 - 07/1989 | -1.969 Dee 
32032 |4 a -9,.134 | 04/2011 - 07/2011 | -2.320 De noto 
32052 (10 o -14.918 | 12/1984 - 09/1985 | -2.416 oO 
32028 16 ds -10.971 | 02/1989 - 07/1989 | -2.173 eos 
sPI-6 [32018 |5 a -10.564 | 03/2011 - 07/2011 | -2.113 o 
32024 |6 e ob -12.945 | 01/2001 - 06/2001 | -2.189 ee /ioso 


23 


32032 |4 eS -8.633 103/2011 - 06/2011 | -2.525 pa 
32052 |15 ió -20.798 | 12/1984 - 02/1986 | -2.437 o oliodo 
32028 |5 bios -10.325 | 03/1989 - 07/1989 | -2.094 eS los 
RDI-6 [32018 16 On -12.613 | 03/2011 - 08/2011 |-2.166 ARO 
32024 (6 a -13.096 | 01/2001 - 06/2001 | -2.243 Desi 
32032 |4 oa -9,.290 | 03/2011 - 06/2011 | -2.622 e 
32052 (6 os -10.257 [06/1974 - 11/1974 | -1.877 PO 
3208 (6 cdo -12.848 | 06/1980 - 11/1980 | -2.141 io 
SPEI-6 |32018 17 al -12.414 | 05/2011 - 11/2011 | -1.881 PS 
32024 (7 o dee -12.949 | 04/1989 - 10/1989 | -1.850 eS 
32032 (9 SN -16.395 | 03/2011 - 11/2011 | -1.822 a 
32052 (40 Deia -65.884 | 06/1985 - 09/1988 | -1.837 Send 
32028 |18 00 -30.437 [07/2011 - 12/2012 |-2.074 O 
spPI-12 |32018  |18 e -38.158 | 02/2011 - 07/2012 |-2.120 e 
32024 |20 a -35.080 | 08/1989 - 03/1991 | -1.754 ba 
32032 23 o -44.154 | 02/2011 - 12/2012 | -1.920 Eon 
32052 |40 pr -61.101 | 06/1985 - 09/1988 | -1.832 dio 
32028 |14 o -24.325 [07/2011 - 08/2012 | -1.835 bo siosa 
RDI-12 132018 (18 pe -36.955 | 02/2011 - 07/2012 | -2.053 aa 
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08/1989 - 08/1989 - 


32024 |20 oca -35.433 | 08/1989 - 03/1991 |-1.772 | ¿37,991 
32032 |23 0 "| -52.654 | 02/2011 - 12/2012 |-2.289 ol i 
32052  |25 e "| -40.618 | 12/1961 - 12/1963 | -1.625 e i 
32028 |12 E "| -22.007 | 06/1980 - 05/1981 | -1.935 iS ; 

SPEI-12 (32018 (18 ol “| -30.817 | 02/2011 - 07/2012 |-1.814 es i 
32024 |24 a "| -38.493 | 06/1989 - 05/1991 | -1.629 a e Ñ 
32032  |26 a “| -50.699 | 11/2010 - 12/2012 | -1.950 o ' 


The results of the SPI and RDI indices were similar at the 3- and 6-month 
scales (SPI-3 and RDI-3; SPI-6, and RDI-6). With SPI-3 and RDI-3 most of 
the drought periods were detected in the months of January to May, while with 
SPEI-3 they were detected from April to July. Regarding SPI-6 and RDI-6, most 
of the dry periods occurred in the months from March to June, while with SPEI- 
6 they occurred from May to September. For these scales, the SPEI detected 
droughts several months after the SPI and RDI, because of the relationship 
that exists between precipitation and temperature during the first months of 
the year. At the 12-month scale, the results of the three indices were similar, 
although in the important periods, in terms of duration, magnitude and 
intensity, SPI-12 and RDI-12 were more similar. According to SPI-12, RDI-12, 
and SPElI-12, major drought periods occurred in 1962, 1969, 1979, 1982, 
1998, 1999, 2000, 2011, and 2012. For most of these years, Campos-Aranda 
(2014) identified droughts at the Fresnillo weather station in central Zacatecas 
using the Probabilistic Precipitation Deficit Index, SPI, and RDI. For those same 
years, Campos-Aranda (2016) also reported annual droughts in a large part of 
Zacatecas using the aridity index. The most important drought period in terms 
of intensity occurred between 2011 and 2012. 


The evolution of the SPI, RDI, and SPEI series at the 3-month scale (short- 
term) had a high temporal frequency of periods of drought and moisture. As 
the time scale increased to 6 and 12 months, the dry and wet periods had a 
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lower temporal frequency and a longer duration (Figure 3). The dry and wet 
seasons were more clearly defined, especially on the 12-month scale. 


3 3 3: 
SPI-3 SPI6 SPr-12 
2 2: 2. 
1 1 1 
50 50 50 
1 -1 1 
2 2 2 
3 -3 3 
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3 3 3 
RDI3 RDI-6 RDI-12 
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Figure 3. Temporal evolution of drought at 3-, 6-, and 12-month scales at 
the El Sauz station (32018). 


Mean drought characteristics 


Table 4 presents a summary of the number of drought events and the mean 
duration of drought periods for short-, medium-, and long-term time scales, 
determined with the three drought indices evaluated. The mean number of 
drought episodes decreased with longer time scales. 
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Table 4. Number and mean duration of the drought events determined by 
SPI, RDI, and SPEI for the period 1961-2012. (The records from this station 
is for the period 1971-2012. Mean duration is shown in parenthesis.) 


Station SPI-3 RDI-3 SPEI-3 | SPI-6 RDI-6 SPEI-6 | SPI-12 el e 
32052 43 44 38 28 30 27 16 14 19 (6.1) 


(2.4) (2.2) (3.2) (3.1) (3.0) (4.1) (6.8) (8.1) 


36 41 47 41 39 34 21 16 
22028 (2.8) (2.7) (2.3) (2.6) (2.7 (2.9) (5.1) (5.8) E 


37 42 45 32 26 34 10 11 
anos (2.9) (2.5) (2.6) (3.0) (3.2) (3.2) (6.2) (5.6) 19:00:80 


42 46 47 38 34 36 13 12 
32054 (2.3) (2.3) (2.1) (3.1) (3.0) (3.1) (8.9) (9.4) 13 (9.1) 


48 49 51 37 36 35 18 19 
GA0HS (2.3) (2.4) (2.2) (2.9) (2.8) (3.0) (5.7) (5.1) Ei 


41 41 48 39 33 33 19 20 
22008 (2.3) (2.3) (2.0) (2.8) (3.0) (3.3) (5.2) (5.0) A 


* 31 34 40 26 26 28 16 19 
32024 (2.8) (2.5) (2.1) (3.1) (3.0) (3.3) (5.4) (4.8) 18 (5.3) 


47 46 44 38 37 39 18 15 
32032 (2.5) (2.6) (2.1) (2.8) (2.9) (2.9) (5.7) (6.7) 13 (8.2) 


44 47 52 34 31 37 19 21 
2030 (2.5) (2.4) (2.0) (3.1) (3.1) (3.1) (4.8) (4.2) Sa 


The average number of drought events at a short-term time scale varied from 
36 - 48 months for SPI-3, from 41 - 49 for RDI-3, and from 38 - 52 months 
for SPEI-3. The mean duration of the drought events ranged from 2.3 to 2.9, 
2.2 to 2.7, and 2.0 to 3.2 months for SPI-3, RDI-3, and SPElI-3, respectively. 
At 70% of the weather stations, the drought event of largest magnitude 
reached the severe drought category with SPI-3 and RDI-3, while with SPEI- 
3, it reached the category of extreme drought at 90% of the weather stations. 
The longest duration of drought periods registered with SPI-3 and RDI-3 
ranged from 5 to 8 months, and from 4 to 7 months for SPEI-3. 


Over the last two decades, the most intense drought events occurred at the 
3-month scale. According to SPI-3 and RDI-3, these events were from 
February to May of 1999 and from January to May of 2011, and according to 
SPEI-3 they occurred from May to July of 2011 and from May to June of 1998. 
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For the 6-month scale, the drought events detected by the three indices were 
relatively similar. Average droughts ranged from 28 to 41, 26 to 39, and 27 to 
39 months for SPI-6, RDI-6, and SPEI-6, respectively. The mean duration of 
these events ranged from 2.6 to 3.1 months for SPI-6, from 2.7 to 3.2 months 
for RDI-6, and from 2.9 to 4.1 months for SPEI-6. With SPI-6 and RDI-6, the 
most intense events reached the category of extreme drought at all the 
stations; with SPEI-6 only 20% reached this category and the rest of the 
stations registered severe drought. According to SPEI-6, the longest durations 
of dry periods were very similar among the stations, ranging from 6 to 9 
months, while according to SPI-6 and RDI-6, the range was greater, from 4 to 
15 months. 


According to SPI-6, RDI-6, and SPEI-6, the northeastern part of Zacatecas had 
major droughts in the 1970s and 1980s, in terms of magnitude and intensity. 
During the 1990s and the first decade of 2000, major droughts occurred in the 
central and southern parts of the state. 


At the 6-month scale, RDI-6 detected the longest drought period (15 months), 
which occurred in the northern region of Zacatecas (station 32052). The 
longest drought detected by SPI-6 and RDI-6 was at station 32001 (northern 
Zacatecas), lasting 13 months, from June 1987 to June 1988. SPEI-6 detected 
a drought for this same period and station, but in two subperiods, one 6 
months long (from July to December 1987) and another lasting 3 months (May 
to July 1988). The longest drought detected with SPEI-6 was 9 months at 
stations 32032 and 32054, from March to November of 2011, in southern and 
west-central Zacatecas, respectively. In this same period, SPI-6 and RDI-6 
detected a period of drought of only 4 months, from March to June of 2011. 


Regarding the 12-month time scale (medium- and long-term), the average 
number of drought events was similar for the three indices, ranging from 10 
to 21 for SPI-12, from 11 to 21 for RDI-12, and from 13 to 21 for SPEI-12. 
The mean duration of the drought events ranged from 4.8 to 8.9 months for 
SPI-12, from 4.2 to 9.4 months for RDI-12, and from 4.6 to 9.1 months for 
SPEI-12. At this scale, with SPI-12 and RDI-12 the most intense events 
reached the category of extreme drought at 70% of the stations, and 30% 
severe drought, while with SPEI-12 only 10% of the stations showed extreme 
drought and the remaining 90% reported severe drought. The longest dry 
periods were similar with SPI-12 and RDI-12, and ranged from 16 to 40 and 
from 14 to 40 months, respectively, while the periods were shorter with SPEI- 
12, ranging from 12 to 27 months. 
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During the 52 years of climatological records analyzed in the state of 
Zacatecas, the SPI detected the most intense droughts on the annual time 
scale, coinciding with the results obtained by Campos-Aranda (2014, 2015), 
who analyzed droughts with SPI and RDI in central Zacatecas. With the 3- and 
6-month time scales, the most intense droughts were detected by SPEI and 
RDI, respectively. 


The most intense droughts at the three time scales occurred between 1998 
and 2011. RDI and SPEI also identified extreme droughts in 1974, 1979, 1982, 
and 1989. Figure 4 shows the spatial distribution of the SPI, RDI, and SPEI for 
the month with the greatest drought at the three time scales analyzed. With 
the 3-month time scale, SPI-3 and RDI-3 identified April 2011 as the month 
with the highest drought intensity. These indices have a spatial behavior that 
is practically equal (Figures 4a and 4b). Most of Zacatecas (82%) suffered 
severe drought and only 8% extreme drought. The most intense drought 
detected by SPEI-3 occurred in June 2011, with extreme drought in 59% of 
the state, mainly in the central part (Figure 4c). At the 6-month scale, April 
2011 was also the month of greatest drought intensity according to SPI and 
RDI; extreme drought occurred in practically the entire state (Figures 4d and 
4e). The SPEI-6 identified July 2011 as the month with the greatest drought; 
most of the territory of Zacatecas underwent extreme drought (49%) and 
severe drought (39%). Extreme drought occurred mainly in central and 
northeastern Zacatecas. At the 12-month scale, September 2011 was 
identified by the three indices as the month with the highest drought intensity. 
SPI-12 and RDI-12 had similar spatial behavior, identifying around 55% of 
Zacatecas with extreme drought (Figures 4g and 4h), mainly in the south and 
west of the state. SPEI-12 identified extreme drought in 33% of Zacatecas, in 
the same regions as did SPI-12 and RDI-12, but in a smaller area. 
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Figure 4. Spatial distribution of the SPI, RDI, and SPEI indices for the month 
with the most intense drought at scales of 3, 6, and 12 months. 


Drought frequency 


Drought frequency was analyzed using the different drought classes (Table 2) 
related to the values of the drought indices (SPI, RDI, SPEI). At the 3-month 
scale, for the drought indices analyzed, the average time (percentage of 
average frequency) of dry and wet periods was 17.1 and 16.9%, respectively. 
Normal precipitation conditions were present 66.0% of the time. SPI-3 and 
RDI-3 showed a higher frequency of moderate droughts (12.5 and 12.2%, 
respectively) than SPEI-3, while the frequency of severe droughts was lower 
for SPI-3 and RDI-3 (4.1 and 4.8%, respectively) than for SPEI-3. Of the three 
indices, SPI-3 and RDI-3 identified a lower number of extreme droughts. 
Figure 5a shows the frequency percentage of dry and wet periods for SPI-3, 
RDI-3, and SPEI-3. With SPI-3 and RDI-3, only the west-central region of the 
study area (station 32054) had extreme droughts (14 and 15 events, 
respectively), while with SPEI-3, all the stations had extreme drought events, 
ranging from 3 to 16 events at stations 32052 and 32032, in northern and 
southern Zacatecas, respectively. 
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Figure 5. Frequency distribution of the values of the SPI, RDI, and SPEI 
indices: a) 3-month scale; b) 6-month scale; and c) 12-month scale for the 
period 1961-2012. 


At the 6-month scale, according to the three multiscale indices, normal 
precipitation conditions were present an average of 67.4% of the time. The 
frequency of dry events was 16.7% and that of wet events was 15.9%. SPEI- 
6 resulted in a higher frequency of moderate droughts (11.0%) than SPI-6 
and RDI-6 (8.7 and 7.7%, respectively). The same occurred with severe 
droughts (Figure 5b). However, SPI-6 and RDI-6 presented a higher frequency 
of extreme droughts than SPEI-6: 3.4 and 3.8% respectively. With SPI-6 and 
RDI-6, 26 and 29 extreme droughts in central Zacatecas were detected 
(stations 32018 and 32030), while with SPEI-6 the highest number of extreme 
droughts detected was 15, which occurred in the north (station 32001). 
Finally, at the 12-month scale, the average frequency percentage of normal, 
wet, and dry events was 66.6, 16.9, and 16.5%, respectively. The frequency 
of droughts with SPI-12, RDI-12, and SPEI-12 was 16.3, 15.8, and 17.4%, 
respectively. The behavior and values of the frequency of the different classes 
of droughts for SPI-12, RDI-12, and SPEI-12 were similar to those detected 
by SPI-6, RDI-6, and SPEI-6. And SPEI-12 presented a higher frequency of 
moderate and severe droughts than SPI-12 and RDI-12. The opposite occurred 
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for extreme droughts, as seen in Figure 5c. With SPI-12 and RDI-12, 23 and 
22 extreme drought events were detected at station 32032 (southern 
Zacatecas), while with SPEI-12, the highest number of extreme drought 
events was 17, and occurred at station 32001 (northern Zacatecas). 


Correlation of the drought indices 


Considering the three time scales and the nine stations analyzed, the SPI and 
RDI were the mostly highly correlated indices, with correlation coefficients (r) 
ranging from 0.920 to 0.996. SPI and SPEI had the lowest correlation 
coefficients, which ranged from 0.336 to 0.982. Figure 6 shows the general 
trend of the relationships between SPI vs RDI, SPI vs SPEI, and RDI vs SPEI 
at the 3-, 6-, and 12-month time scales, with dispersion diagrams for the 
series of indices at El Sauz station (32018), located in the central region of 
Zacatecas. 
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Figure 6. Dispersion diagrams for the SPI vs RDI, SPI vs SPEI, and RDI vs 
SPEI series for the El Sauz station (32018), period 1961-2012. 


At the 3-month scale, the correlation coefficients (r) for SPI-3 vs RDI-3 ranged 
from 0.920 to 0.969, while SPI-3 vs SPEI-3 and RDI-3 vs SPEI-3 ranged from 
0.336 to 0.804 and from 0.583 to 0.867, respectively. In general, the 
correlations for the 6-month scale were slightly higher than for the 3-month 
scale, ranging from 0.943 to 0.964 for SPI-6 vs RDI-6, from 0.410 to 0.834 
for SPI-6 vs SPEI-6, and from 0.643 to 0.906 for RDI-6 vs SPEI-6. At the 12- 
month scale, the correlation coefficients ranged from 0.952 to 0.996 for SPI- 
12 vs RDI-12, from 0.926 to 0.982 for SPI-12 vs SPEI-12, and from 0.952 to 
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0.993 for RDI-12 vs SPEI-12. For the scales and stations analyzed, the smaller 
the scale, the lower the correlation between indices (Table 5). At the 3- and 
6-month scales, the lowest correlations were in the north of the state (station 
32052) and the highest in the south (station 32032). 


The Spearman correlation coefficients (Table 5) were similar to those of the 
Pearson correlation coefficients. As for coefficient r, at the 3- and 6-month 
scales the SPI vs RDI indices had higher values of r,, while SPI vs SPEI had 
lower values. For the 12-month scale, the RDI-12 vs SPEl-12 and the SPI-12 
vs RDI-12 had the highest r, coefficients. The similarity between the r and 
r, coefficients increased as the scale increased. 


Table 5. Correlation coefficients r and r, for the SPI, RDI, and SPEI drought 


indices. 
Statio 
n 32001 32003 32028 32018 32032 32054 32052 32030 32024 
Coefficie 
Index nt r Ys r Ys r Ys r Ys r Ys r Ys r Ys r Ys r Ys 
SPI-3 RDI-3 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 
33 28 55 48 48 50 47 46 69 69 49 48 20 26 34 40 59 56 
SPEL-3 0.4 0.3 0.6 0.6 0.5 0.5 0.6 0.6 0.8 0.8 0.7 0.7 0.3 0.2 0.7 0.7 0.6 0.6 
65 82 71 75 62 23 36 33 04 26 65 75 36 41 66 73 33 22 
RDI-3 | SPEI-3 0.6 0.6 0.7 0.8 0.7 0.7 0.7 0.8 0.8 0.9 0.8 0.9 0.5 0.5 0.8 0.8 0.7 0.7 
84 54 93 34 25 17 78 11 67 10 62 06 83 31 62 99 66 83 
SPI-6 RDI-6 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 
43 33 58 31 47 35 49 29 62 39 52 36 44 41 43 11 64 41 
SPEL-6 0.4 0.4 0.6 0.7 0.5 0.5 0.6 0.6 0.8 0.8 0.8 0.8 0.4 0.3 0.7 0.8 0.6 0.6 
99 45 83 05 60 68 37 55 34 75 01 21 10 31 71 00 54 77 
RDI-6 | SPEI-6 0.7 0.7 0.8 0.8 0.7 0.7 0.8 0.8 0.9 0.9 0.9 0.9 0.6 0.6 0.8 0.9 0.7 0.8 
16 04 12 91 50 82 03 58 06 70 01 45 43 00 86 48 90 60 
SPI- RDI-12 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 
12 96 96 96 96 87 88 96 95 52 36 91 90 96 96 80 77 91 89 
SPEI-12 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 
64 80 82 86 50 51 81 85 26 12 76 81 27 40 52 59 62 60 
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RDI- 
12 


SPEI-12 


0.9 0.9 | 0.9 0.9 0.9 | 0.9 0.9 0.9 [0.9 0.9 0.9 | 0.9 0.9 0.9 | 0.9 0.9 
78 92 93 97 86 85 93 97 92 97 91 98 52 60 91 97 


0.9 
87 


0.9 
89 


In this study, SPI and RDI were more suitable for characterizing meteorological 
drought. According to Gocic and Trajkovic (2014) and Xu et al. (2015), RDI 
and SPI are more suitable indices than SPEI for characterizing drought in arid 
and semiarid regions. Our results coincide. 


Moreover, selection of the probability distribution function for the analysis of 
drought indices may affect the results, underestimating or overestimating 
drought events. Therefore, it is helpful to use a frequency analysis to 
determine the distribution function that best fits the precipitation, 
precipitation/evapotranspiration, and precipitation — evapotranspiration series 
when evaluating the SPI, RDI, and SPEI indices, respectively. 


Conclusions 


The multiscale meteorological drought indices used in this study made it 
possible to characterize the droughts that occurred over the last 50 years in 
the state of Zacatecas, Mexico, providing the number, beginning, end, 
duration, magnitude, intensity, and category of the events. 


The three drought indices consistently detected drought events, particularly 
those of greatest magnitude, regardless of the time scale of the index. They 
also enabled conducting a complementary analysis of the spatial and temporal 
variability of the droughts. 


The events detected with the three indices had a lower frequency and longer 
duration as the time scale of the analysis increased. The three scales presented 
drought events characterized as moderate and extreme, which were of short 
and medium duration at the 3- and 6-month scales, and of medium to long 
duration for the 12-month scale. The RDI detected a higher number of 
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droughts at the 3-month scale, while SPI and SPEI did so at the 6- and 12- 
month scales, respectively. 


The spatial analysis of the droughts shows that the areas that were most 
affected by this phenomenon were the northwestern, central, and southern 
parts of the state. The principal rainfed and irrigated agricultural regions of 
Zacatecas are located in these areas. 


The Pearson and Spearman correlation coefficients were higher as the scale of 
the indices increased. The r and r; coefficients confirmed that there was 
greater similarity between the results of SPI and RDI for the 3- and 6-month 
scales. At the 12-month scale, the performance of the three indices was 
similar, and the results of RDI and SPEI had greater similarity. Therefore, if 
information on temperature is available, it is recommendable that the RDI be 
used to analyze droughts in semiarid regions, since it incorporates other 
climatological variables that have an important role during the development of 
a drought event, such as evapotranspiration. Nevertheless, in absence of 
temperature information or for ease of calculation, for the scales of analysis, 
SPI is a good option for estimating drought characteristics. 


Characterization of historical droughts is important for generating information 
related to water deficits during periods of drought, since it can contribute to 
designing policies aimed at reducing damage from future drought events. 
Given the importance of the agricultural sector in the state of Zacatecas, the 
findings of this study may support decision-making mainly regarding crop 
patterns and optimal water allocation in Irrigation District 034 in the state of 
Zacatecas. 
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